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Abstract. We solve the dynamics of large spherical Minority Games (MG) 
in the presence of non-negligible time dependent external contributions to the 
overall market bid. The latter represent the actions of market regulators, or other 
major natural or political events that impact on the market. In contrast to non- 
spherical MGs, the spherical formulation allows one to derive closed dynamical 
order parameter equations in explicit form and work out the market's response 
to such events fully analytically. We focus on a comparison between the response 
to stationary versus oscillating market interventions, and reveal profound and 
partially unexpected differences in terms of transition lines and the volatility. 



PACS numbers: 02.50.Le, 87.23.Ge, 05.70.Ln, 64.60.Ht 
1. Introduction 

Minority Games (MG) jHll] are simple mathematical models designed to elucidate and 
explain the origin of the nontrivial macroscopic fluctuation phenomenology observed 
in real markets, on the basis of so-called inductive decision making by large numbers 
of interacting agents [Sj. Their great strength is that they can be solved analytically 
using methods from the statistical mechanics of disordered systems, in particular with 
generating functional analysis techniques [11 O [6j . We refer to the recent textbooks 
|7j ^8J for historical backgrounds, the connection between MGs and real markets, 
details on mathematical methods, and full references. Now that the standard MGs 
have been solved and understood satisfactorily, attention must turn to generalizing 
the mathematical technology developed so as to apply to models that are more 
realistic economically. In this paper we try to contribute to this aim by studying 
analytically the dynamical response of MG markets to non-negligible events in the 
outside world that impact on the overall market bid (which in MGs is a proxy for 
the asset price). These events could be accidental (e.g. natural disasters, changes in 
resource availability), political (e.g. election results, major management incompetence 
or corruption scandals), or interventions by market regulators. Such ingredients are 
incorporated easily into the fabric of MG-type models, by simply adding to the 
overall market bid time-dependent external terms; the tricky stage is to work out 
mathematically the consequences of such terms. The simplest case is that where 
the external bid term is stationary. Here only minor modifications of the standard 
formalism are required; see e.g. [S] (using the replica method) or [lOj (using generating 
functional analysis). For non-stationary external bid contributions, in contrast, one 



Market response to external events and interventions in spherical minority games 2 



is generally forced to either resort to numerical simulations (see |10|). or turn to 
those generating functional analysis versions that involve explicit representations of 
the overall bid process, as developed for MGs with real market histories [UlIIl]. In 
the latter studies only infinitesimal external bid perturbations have been considered 
so far, in view of the complexity of the formalism. There is, however, one special 
class of MG versions where adding time-dependent external bids does not lead to 
serious mathematical complications: the spherical MGs as introduced in [13l [14]. 
We show in this paper that here it is still possible to derive fully explicit and exact 
equations for time dependent order parameters (correlation and response functions); 
time-dependent external bids are found to be mathematically harmless, and one can 
even allow for (partial) market impact correction as in [15j . We focus on comparing 
the long-time solution of these equations for stationary external bids to those found 
for oscillating external bids. As expected, these two cases are found to generate very 
different macroscopic consequences, in terms of phase diagrams and in terms of the 
volatility. Some of these are intuitively clear, such as the profound impact of oscillating 
market interventions on the volatility, but some are not at all intuitive, such as the 
independence of the phase diagram in the case of oscillating external bids on the 
amplitude of the intervention (in contrast to stationary bids). All our theoretical 
results are tested against numerical simulations, and find excellent confirmation. 



2. Definitions 



In MGs one considers N agents, labeled usually by Roman indices i = 1, . . . , N. At 
each time step t S {0, 1,2,.. .} of the game each agent i submits a trading action, a 
'bid' bi(t), in response to public information fj,{t) which in fake history MG versions 
is chosen randomly and independently from the set {1, . . . ,p} (where p = aN with a 
finite). The rescaled total market bid in the game at time t is subsequently defined as 

A(t)=A,(t) + -^^6,(t) (1) 



Here Ae{t) could represent e.g. random market perturbations, actions by market 
regulators, or other external events (natural, social, political, etc.) that can change 
the overall asset demand in the market directly. Each agent i has S 'look-up table' 
strategies R*"" = (i?^°, . . . , R^) e {-1, 1}^, with a = 1, . . . , 5. If agent i decides to use 
strategy a at time t in the game, his bid at that stage will be bi{t) — R^^^^ty In the MG 
an agent i finds himself winning at time t if his decision hi{t) turns out to be opposite 
in sign to the total bid, i.e. if hi{t)A{t) < 0. All agents monitor the performance of 
their strategies, in order to decide which one to use. To this end, they assign points 
to each of their strategies based on the update rule 

P^a{t + 1) = P^ait) V^^) [A{t) ^ [iT^^f - R^^^)] (2) 



N 

where is a learning rate (which sets the unit time scale), and where the first 
minus sign implements reward for minority decisions. The term proportional to 
K S [0, 1] represents a (partial) correction by individual agents of A{t) for their 
own contribution, as in [15 . The strategy played by agent i at time t is then 
ai(t) — a,TgT[ieLXa£{i,...,s} Piait). In this paper we limit ourselves to S* = 2; viz. two 
strategies per agent. It is now sufficient for agents to keep track only of the differences 
gj(t) — ^[pii{t) —pi2(t)]- Upon also replacing in the right-hand side of ^ the random 
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at each step t by an average over all possible values (the so-called batch version 
of the game), viz. RHS[^(i)] p^^ J2fi<p R.HS[/z], and with r/ = 2a\/lvQ one finds 

q.{t + 1) = q,it) + e.it) - ;4 E (^''(*) - ^^'^'(^^^■O 

A'<P 

(t) = A, (t) + + ^ </>^. (t)^; (4) 

with^'^ = ^{R^-R^),u;^ = =iV-V2^^^,f,and 0,(t) = sgnfe(0]. 

We have also added a perturbation field 9i{t) to define response functions later. Finally 
we introduce a spherical constraint into the model, of the type proposed in [Tl] : we 
replace the previous (so-called batch MG) definition (j)i{t) = sgn[qi{t)] by 

Mt) = q^{t)/X{t), X{t) - ^ q^itf]^/^ (5) 

i 

The above spherical MG version is a generalization of f4j. It is not unique; the 
alternative spherical MG in I3J differs from the present formulation in at what stage 
and for which variables the relevant nonlinearities are replaced by pseudo- linear laws. 



3. Generating functional analysis 

3.1. Derivation of exact order parameter equations 

The generating functional analysis (GFA) method is based on the calculation of the 
generator Z[ip] = (exp[— i^-^ Vi(i)0i(O])j by interchanging the averages (...) over 
paths (which here, in the absence of decision noise, reduce to averaging over initial 
conditions) and over the disorder (i.e. the strategies). It has become the standard 
tool to study MG dynamics, so we may refer to e.g. [H] for technical details. For 
N ^ oo and upon choosing 9i{t) = 9{t) for all i, the method leads to the following 
self-consistent equations for two-time correlation- and response functions Cw and Gw , 
written in terms of averages (. . .)* over an effective single agent process: 

t^t': Cw = mm')),, Gtf =dm))*/d0{t') (6) 

A further order parameter X{t) is to be solved from Ctt — 1, reflecting the spherical 
constraint, and causality ensures that Gu' = for all t < t' . The effective single agent 
process is defined by the following stochastic equation, with (f){t) = q{t)/X{t): 

q{t + 1) = q{t) + e{t) - a E [ (I+G)" - + V^v{t) (7) 

t'<t 

Here ri{t) is a zero-average Gaussian noise, characterized by {ri{t)ri{t' )) = Y.[Ag]tt', 

EKlu' = [(I + G)-'D[A,]iI + G^)-%' (8) 
D[A,]tt' = 1 + Gtt' + 2A,{t)A,{t') (9) 

The effective Gaussian noise replaces the statistics of the original N agents by evolution 
uncertainty for one effective agent. We can now appreciate the advantages of the 
spherical version. The effective equation ([7]) is linear in (j){t), which allows us to derive 
fully explicit dynamical equations for the order parameters. The only nontrivial step 

t The Vn ensures that the relevant time scales are 0{N^); the factor 2a leads to simple equations. 
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in this derivation is to show via integration by parts that {ri{t)<f){t'))i, = (S[yle]G^)tt/. 
For 0{t) ^ Q the final result is 

\{t+l)Ct+i,v - [m + HCtt' = + G)-^D[A,]{1 + G^)-^G%, 

~a[{\ + G)-^C]u' (10) 

\{t + \)Gt+i,v - [A(0 + aK]Gu' = <5tt, - a [(I + G)-^G] (11) 

By solving these equations for the kernels {C, G}, together with the condition Ctt — 1 
from which to extract A(i), we can explore the dynamics of the original MG model for 
N — > oo, since the physical meaning of {C, G} is 



^E(^^W<^«(^')> (12) 



^"'= J™ ^jY.^mt))ido.{t') (13) 

i 

The advantages of the spherical MG are that one can derive an exact formula for the 
volatility (see below), and that the explicit nature of its order parameter equations 
allows us to analyze the effects of the bid perturbations {t) much more effectively 
that in ordinary MGs. 

3. 2. Bid average and fluctuations 

The statistics of the overall bids (t) in (|4]) can once more be extracted from a suitable 
generating functional, namely Z[ip] — (exp[iV2^^j (^^(t)A^(t)]). It generates the 
relevant moments of the overall bids via differentiation, e.g. 



dZ[cp] 



(A^'(t)) = ^ lim ^ (14) 

^ ^ " V2 v^O d^^{t) ^ ' 



1 d'Zy] 



lAt^it^A^it')) = - - lim - — , , , (15) 

Following the familiar steps of the generating functional analysis technique leads us 
back to the previous saddle-point problem, but now we obtain an additional expression 
for the bid moments. We refer to [SI [10] for full details of such calculations, and limit 
ourselves here to giving the final result: 

WW) = ^(I + G)-,Me(0 (16) 

t' 

WifWW = WW WW + ^'5m-[(I + G)-'Do[A,]{\ + Gt)-i],,, (17) 

where Do[Af.]tt' — 1 + Gtt' (note that Do[^e] still depends on Af.{t) via the kernel C). 
Clearly, in the presence of finite external bid perturbations the system is no longer 
guaranteed to evolve towards a state with zero-average bid statistics. We can now 
define a fluctuation volatility by the following expression: 



al^ hm ±^^|([A^(t)p)_(AM(t)) n 

r—*oo TV ^ — ^ ^ — ^ 1^ J 

^ t< 



= hm ^^[(I + G)-ii^oK](I + Gt)-i],, (18) 

t<T 
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Note, however, that in the presence of non-stationary Ai,{t), even in the absence of 
anomalous response and upon assuming self-averaging with respect to the disorder 
it will generally no longer be true that aa is identical to the disorder-averaged 
conventional volatility a as defined by 

2 



lim — VVM^mi^- [ lim — VVA^(t) 

t<T H t<T fJ. 



(19) 



The difference between a and aa reflects bid oscillations which are deterministic and 
therefore excluded from (Tfl. We will derive an exact relation between the two later. 



4. Time translation invariance with constant or oscillating external bids 

In the remainder of this paper we focus on the choices Ae{t) — A and Ae{t) = 1)*, 
where the asymptotic consequences of bid perturbation are most easily quantified. 
Experience with previous MG versions suggests that there will be two types of players 
in the stationary state: 'frozen' agents, with qi(t) growing linearly with time, and 
'fickle' agents, where qi{t) does not diverge with time. In view of this we consider two 
types of solutions with respect to the spherical constraint parameter X(t): a regime 
where limt^oo X{t) — X (finite), and another regime where X{t) ^ oo a.s t ^ oo. 



4.I. General formulae 

For Ae{t) — A(— 1)'»* with C £ {0, 1} our equations (jlOllip will have time-translation 
invariant (TTI) solutions, since here jS]) gives D[Ae]tt' = 1 + Cw + 2i2(_i)C(t-t'). 
Upon assuming Cu' ~ C{t — t') and Gu' — G{t — t'), so that the same is true for 
D[Ae] and E[yle], our equations pOllip then reduce for finite values of t to 

G{t+1) - ^boGit) ^ a^i{^[{l+ G)-^D[A,]{1+ G^)-^G^]it) - [(1+ G)"' q(t)| 

(20) 

G{t+1) - 7AoG(i) = Vijjto - a [(1+ G)-'G] (t)} (21) 

with 

^/>o = lim [A(t-l) + aK]/A(i), Vi = 1™ 1/A(t) (22) 

t — >CXD t — >00 

If X{t) ^ A for t ^ oo one has ipo — {X + aK)/X and ipi — 1/A, whereas if X{t) — > oo 
one has ipo = 1 and tpi = (note: since strategy valuations can diverge at most 
linearly with time, the same must be true for X{t) so limt_oo A(i — 1)/X{t) < 1 is 
ruled out). There is, however, an important subtlety. It is not clear that in the case of 
diverging X{t) one can use equations (|20I21[) to calculate observables such as x or x, as 
this requires that the limits t ^ oo in x = X]t<T i — > oo in (|22p commute. If 

they do not commute (as will be the case), this means that the time it takes to evolve 
from initialization to a state with {ipo,tpi) — (IjO) diverges, so such states will in 
practice not be observed. Equations (|20l2ip invite us to switch to Fourier transforms, 

at) = r ^e'-'Cic), C{uj) = e--*C(i) (23) 

The Fourier transform of the kernel ([5]), in particluar, is seen to be 

i)[Ae](w) = 27r(5(w) C{u:) + 47ri2j(t^ _ (7^) (24) 
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One expects the relevant static observables to include the integrated responses to 
static and oscillatory fields, viz. x = J2t = ^(^) ^^"^ ^ — Z]t(~l)*^(^) = G(7r). 
With these definitions we can write (|20l2ip as 

Ciuj) [|l + G(c^)|2(e'" - Vo) + + G(-to) - G{uj)] 

^2TTaTpiGiLu)[5{Lo) + 2A'^6{Lu-CTT)] (25) 

(e'"" - ^o)Giuj) = V^i [l - aG(cj)/[l + G(cc;)]] (26) 

We must also enforce the spherical constraint G(0) ~ 1, which gives \AijjC[lo) — 2-k. 
The right-hand side of (I25p dictates that the solution must me of the following form, 
where we have built in the spherical constraint, with cq G [0, 1]: 

C{lj) = 27r[co(5(w) + (1 - cq)5{uj - tt)] (27) 

Equations (|25I26[) then reduce to the following set: 

co[a^i + (l+x)'(l-V^o)] = a^ix[l + 2A^5^q\ (28) 

(l-co)[aV;i - (l+x)'(l + V^o)] = 2a7^ixi'^ci (29) 

(i-Vo)x(i + x) -V^i(i + x-ax) (30) 
-(i + V^o)x(i + x)-V'i(i + x-«x) (31) 

We can next work out the bid statistics in TTI states. Here formulas (|16ll8p are 
seen to give the following, which shows (as expected) that the effects of stationary or 
oscillatory external bid perturbations are generally reflected in the bid averages: 

\^-YMm-^. hm i^(-l)*(A^.^ (32) 

4-2. Solution with diverging constraining force: fully frozen states 

In states where X{t) ^ cx3 for i — > oo one has ('0o,V'i) = (l^O), so equations (|20l2ip 
reduce to C{t) ~ 1 and G{t) — for all t, indicating a fully frozen microscopic 
state with cq — 1. The simplest potential solution of our remaining equations is 
that where we subsequently assume the time limits in the definition of {x, x} to 
commute with those in (gH). It then follows from (|28I29I30I31|) that ^ = X = 0, 
and the volatility ([55)1 would become — 1. Unfortunately, the relevant limits do 
not commute, which follows from a more careful analysis of the effective single agent 
process. Although still cq = 1, the calculation of x ^nd x has to be done explicitly. 
Upon defining Q = liiat^oo q{t)/t, A = limi^oo A(t)/t, rj = limt^oct^^J2t'<t'^i^')' 
and 9 = linit^oo t^^ X]t'<t and upon writing ^ in integrated form, one finds 

~9 + V^V (34) 



^ l + (a/A)[(l+x)-i-^] 

This gives Qirj), where ry is a zero- average Gaussian random variable with (?7^)* = 
lim^^o, t-i Ess'kM^Ms')) = il + co + 2A\o)/{l + x?_ = 2(1 + AH^^)I{1 + x? ■ 
The susceptibility now follows from x = ^^^d{Q{T]))*/d9 = [K^/a\^^ {dQijj)/ drf)i,, 
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and A follows from the spherical constraint linit^oo(0^(^))* = 1; which translates into 
(Q'^iv))* = A^. So we find the following two coupled equations, respectively, 

^ = A[l + (a/A)[(l + x)-i-«]] ^^^^ 



^ ^ (l + x)[l + (a/A)[(l + x)-i-/«]] ^^^^ 
To calculate x we require an infinitesimal oscillating field 9{t) — (— 1)*6', and use 

X= hm^ lim i5](-l)*(||l). = (37) 

by virtue of q{t)/X{t) — > Q{fi)/K as t — > cx). Solving the above equations for x a-nd A 
and substituting the result into (|33|1 gives 

X - , = (38) 

/2a(l + i25^Q)_l 



A = -l-a(l-K) + x/^(3 + 2yl2,5^o)/V2(l+^^'5co) (39) 



2 



(l-l/V2a(l+A2<5co)j (40) 

Ergodicity breaks down when x ^ oo, i.e. at ac,i{A) = i(l + A^Sqq)'^ . Secondly, 
we see that the above solution breaks down when \{t) no longer diverges with time, 
which happens when [1 + a(l - k)]^ = \a{?, + 2A^5i^oYl{l+A^5Qa), i-e. at 



~ R-2{l-n)±^R^-^l-K)R {i + 2AH^^f 

ac±[A,K) = — , R= (41) 

' 2(l-/t)2 2 + 2v42^C0 

The solution Ofc- is unphysical as it obeys ac- < ac,i, which implies that it occurs 

in the nonergodic regime where the above formulae are no longer valid. The solution 

Qc,+ is relevant as it obeys ac.+ > olc,i- Thus the regime of diverging \{t) is bounded 

on the left (as a function of a) by a non-ergodicity transition at a = ac^i{A) and on 

the right by a transition to states with finite \{t) at a = ac,+ {A, k). The transition 

value ac,+ {A, k) increases monotonically with increasing k; it is minimal at k = 0, 

taking the value ac,+ {A,Q) = ^[R—2 + ^/R"^ — AR\ , and diverges at k = 1. 

4-3. Solutions with finite constraining force 

In the alternative scenario limt_^oo ^{t) = A e IR we substitute into (|28l29l30l3ip the 
values tpQ = (A + aK)/A and V'l = If we also write A as A = ^a{-^~ n), to simplify 
our equations further, we find 

X + 2xAH^^ 2xAHa 
1-k{1 + xY l-7(l+x)2 
ai^X^ + X[l + a(K - 1)] + 1 = (43) 
a7x' + X[l + a(7 - 1)] + 1 = (44) 
Solving the last two equations for x x gives two possibilities for each: 

a(l-K)-l± J\a{l-K)-lV - Aan , , 

X± = 7y 45) 



q(1-7)-1± V[a(l-7)-l]"-4^ 
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It follows from the demand x G IR that solutions with a finite constraining force exist 
only if [q;(1 — k) — 1]2 — AaK > 0, that is, only if 

Furthermore, upon investigating the limit k ^ (no impact correction) in (|45p we 
find limK^o X- = ^/i'^^^) ^-nd limK^o X+ = oo? so the physical solution must be x- ■ 

which is finite for all k > 0, and positive for a>I/(l — k). To find the static order 
parameter cq we need only combine (|48p with the first equation of (j42p . The remaining 
two equations, viz. (j46[) and the second equation of (|48p . serve only to determine the 
quantities x and 7. Substituting ([48]) into (|42|) gives 



^ _ 2a(l + 2A2j^o)[a(l-«)~l - V[a(l-^)-l]2~4^] 

^" 4a2K - [a(l + K)-l - ^[a(l-K)-l]2 - Aan]^ 

We extract from this that cq — — a — a2(l+2K) +0{a^) for small a, and that cq i —00 
for a t (5c,-- Hence the solution (|48|49p is unphysical in the regime a G [0, cic-]. This 
leaves the regime a > ac.+. Here we are sure that x > (since 1/(1 — k) < ac,+), 
and we find cq = l/a(l — k) + 2K/a'^{l — k)'^ + 0{a~^) for a — > 00, and that cq 1 00 
for a i ac.+ - Hence the solution (|48l49p is physical in the regime a S [fie, 2, 00), where 
'5c,2(^e,K) > 1/ {1~ ^/k)'^ is defined by the condition cq = 1. This last critical line is 
most easily derived indirectly, by first using (j42p to write it as a condition on x, giving 



'(1 + 2^25^0)' + 8k(1 + A^S^o) " 1 - 2^25^0 

Xc,2 = ^ 1 (50) 

Zk 

This is then inserted into (|43)) and the result is easily solved for a, giving 

ac2(Ae,K) = = ^ — = (51) 

^ ^ ^ (S- 1- 2^2,5^0 -2«;) (3 + 2^2,5^0 -S) 

{1 + 2A^S^o)^ + 8k{1+A^S^o) (52) 

After (tedious) reworking of this result, using the definitions of S and R and the 
identity S = [2(1 + A'^S(^a){R — 4(1 — k))Y^^ to connect the various short-hands, one 
finds that this critical line ([5T|) marking cq = 1 with finite X{t) is identical to the 
previous critical line (|^T|l . For A = we recover the corresponding formula of [H] : 



ac,2(0,K) = — —2 (53) 



4k + 5 + 3VTT8k 
4(1 -k)2 

Let us finally turn to the susceptibility x, which we need to calculate the volatility 
(|33p . and which is to be solved together with 7 from (|44| and the second equation in 
(|42p . One finds upon eliminating 7 that 

X± = ^= (54) 



l±^a[l + 2A25^i/(l-co)] 

Since the present solution regime always has a > 1, we see that the two distinct 
solutions x± represent in-phase and out-of-phase responses of the market to the 
oscillating external bid input, and that neither will be able to diverge. The volatility 
(|33)) now follows upon substituting the result ([Si]) , together with (|48)) and (|49|) . 
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5. Phase diagrams 

We can now summarize the picture obtained by analyzing the time-translation invari- 
ant states in terms of a phase diagram, where the control parameters are a, k, and A. 
The three phases of the system are the following: 

a € (ac,2,oo) : oscillating phase (O), cq < 1, x < oo, x 7^ 0, A(t) < oo 
a € (q!c,i, Q!c,2) : frozen phase (F), cq — 1, x < oo, x = 0, X{t) -^oo 

a € [0, ac,i) : anomalous phase (A), cq = 1, x — oo 

The precise dependence of the transition lines ac,i and ac,2 on the two remaining 
control parameters {A, k) is furthermore controlled strongly by whether the external 
bid perturbation is static (C = 0) or oscillatory (C = 1). In the former case one finds 

o^cAA) =i(l + i2)-i (55) 



R{A)-2{l-K)+dR\A)-A{l-K)R{A) {i + 2a^Y , ^ 

Here the frozen (F) phase will consistently grow with increasing values of A, due to 
both its left boundary ac^A) decreasing with A and its right boundary ac.2{A,K) 
increasing with A. This is easily understood: a static external bid perturbation 
diminishes the impact of the internal bid contributed by the agents, and makes it 
easier for the agents to find a suitable strategy that will land them in the minority 
group. At K = (no self-impact correction) the transition point Q!c,2 simplifies to 
ac,2(^,0) = 2(1 + A^) = l/ac^A). For oscillatory perturbations Ae{t) = ^(-1)*, in 
contrast, the situation is very different: here the relevant transition lines are strictly 
independent of the amplitude A, and are obtained by evaluating the above formulae 
at ac,i(0) and ac,2(0, k), viz. 



1 , , 5 + 4k -f 3^1 -t- 8k 
ac,i = 2^ ac,2(K) = 4^(1 -k)^ (^'^^ 

However, although in the case of oscillatory perturbation the phase diagram is 
independent of A, in the oscillating (O) phase there will be a significant dependence 
on A of the volatility. In all cases we find that the effect of self-impact correction is 
to strengthen the frozen (F) phase, with the oscillating phase vanishing altogether for 
full impact correction (i.e. for k = 1). 

It follows from expressions ([5^ that throughout the phases F and O the system is 
not able to compensate fully for stationary external bid perturbations; since x is finite, 
the average bid will be nonzero for any nonzero perturbation amplitude A. Only at the 
transition fine ac,i{A), where nonergodicity sets in and x diverges, does the bid average 
vanish. This is the situation that was also encountered in ordinary MGs [101 ■ In our 
present spherical MG, however, we are able to inspect also the impact of oscillating 
external bid perturbations, and find that with such perturbations the situation is 
worse: since x remains finite everywhere, we will always retain an oscillating (and 
predictable) nonzero overall bid average, even if we enter the nonergodic regime. 
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Figure 1. Phase diagrams for static external bid perturbations, Ag(t) = A, 
exhibiting an oscillating phase (O), a frozen phase (F), and an anomalous phase 
(A). The latter is characterized by x = oo- Dashed lines: the F— +A transition 
ac,i- Solid lines: the F^O transition ac,2- Left diagram: transitions in the (a, k) 

plane. Here wc show the lines for A = 0,1, 2, 3, 4 (left to right in the case of F— >0, 
right to left in the case of F^A). Right diagram: transitions in the {a. A) plane 
for K = (no self-impact correction), where cvc, 2(^,0) = l/ac,i{A). For static 
external bid perturbation, increasing A always increases the F phase. 



A 




1 

a 



Figure 2. Phase diagrams for oscillating external bid perturbations, Ae{t) = 
A(— 1)', with the same definitions of phases and transitions as in the previous 
figures. For oscillating bid perturbations the phase diagram is independent of the 
perturbation amplitude A. 
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K = K = 0.25 K = 0.5 




a a a 



Figure 3. Macroscopic observables for static external bid contributions Ae{t) = 
A. Top row: frozen correlations cq for /t = (left), k = 0.25 (middle) and k = 0.5 
(right); with in each picture the values as measured for A S {0,1,2,3,4,5}. 
Bottom row: volatility a for k = (left), k = 0.25 (middle) and k = 0.5 
(right); with in each picture the values as measured for A £ {0, 1, 2, 3, 4, 5}. In all 
figures the connected markers indicate simulation results (full circles: biased initial 
conditions; open circles: unbiased initial conditions), with solid lines showing 
the corresponding theoretical predictions. The vertical dashed lines mark the 
transition values ac^i{A) (left, transition from the nonergodic regime with finite 
spherical constraint force regime for small a to the frozen regime) and «c,2(^, i^) 
(right, transition from the frozen regime with infinite spherical constraint force to 
the large a ergodic regime with finite spherical constraint force). 



6. Tests against numerical simulations 

In this section we test our predictions regarding the values of the long-time order 
parameters and the locations of phase transition lines, for stationary and oscillating 
external bid contributions, respectively, against numerical simulations. All simulations 
were carried out with systems of size N — 3000, for both unbiased {qi{Q) — ±10"'') 
and biased {qi{0) = ±1) random initializations. Observables were always measured 
over a duration of 2000 batch iterations, following an equilibration stage of 1000 batch 
iterations. When comparing theoretical predictions to simulation measurements it 
became immediately clear that of the two possible saddle-points in ((54|) the physical 
one is x+: the in-phase solution with high volatility; henceforth all theoretical 
predictions will correspond to this saddle-point. 
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Figure 4. Macroscopic observables for oscillating external bid contributions 
Ae{t) = 1)*. Top row: frozen correlations cq for k = (left), re = 0.25 
(middle) and n = 0.5 (right); with in each picture the values as measured for 
A e {0,1,2,3,4,5}. Bottom row: volatility ct for re = (left), re = 0.25 
(middle) and re = 0.5 (right); with in each picture the values as measured for 
A £ {0, 1, 2, 3, 4, 5}. Note the different vertical volatility scales compared to the 
previous figure. In all figures the connected markers indicate simulation results 
(full circles: biased initial conditions; open circles: unbiased initial conditions), 
with thick solid lines showing the corresponding theoretical predictions. The 
vertical dashed lines mark the transition values a^^i = ^ (left, transition from 
the nonergodic regime with finite spherical constraint force regime for small a 
to the frozen regime) and Qc,2{k-) (right, transition from the frozen regime with 
infinite spherical constraint force to the large a ergodic regime with finite spherical 
constraint force). 



The volatility a that is measured in simulations is the conventional one given 
in pop , for which we still have to derive an analytical expression. Upon combining 
definition (fT5|) of an, for which we have derived prediction with (following 
any given initial state, so that the brackets (. . .) are irrelevant, and assuming self- 
averaging of both volatilities over the disorder realization for N ^ oo) we see that 

^ t<T II ^ tt'<T IJ.f 

t<T t' 
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Ae{t) = A A,{t) = A{-iy 




0.01 0.1 1 10 0.01 0.1 1 10 



a a 

Figure 5. Examples of measured values (connected markers) versus predicted 
values (thick solid lines) of the spherical constraint force A, as measured for 
A £ {0,1,2,3,4,5} with k = 0. The vertical dashed lines mark the transition 
values ac,l (left, transition from the nonergodic regime with finite spherical 
constraint force regime for small a to the frozen regime) and etc, 2 (right, transition 
from the frozen regime with infinite spherical constraint force to the large a 
ergodic regime with finite spherical constraint force). The force A is is indeed 
seen to diverge at the predicted values of a. 



-[^J2J2i^ + G)u'Aeit')]'} (58) 

t<T t' 

In time-translation invariant stationary states with Ae{t) E {A, A{—lY} this becomes 
a'^ai + J2il + G)-\s){l + Gr\s') 

ss' 

X ^lirn {-^Ae(i-s)A,(t-s') - ^ M^- s)Ae{t' ^ s')^ 

_ 2 , f iV(l + x)' if A,{t)=Ai-lY 

+ 1 if A,{t)^A ^^^^ 

The results are shown in figures[3llll and[5l One observes excellent agreement between 
theory and numerical experiments, both in terms of the values of the observables and 
the locations of the O— s-F and F— >A transitions. In particular, we find confirmation 
of the somewhat surprising prediction that, although in terms of market stability 
the oscillating external bid term is a disaster for the market (the market always 
responds in-phase to the oscillation, i.e. making the deviations from market efficiency 
consistently worse) , the locations of the phase transition lines remain always identical 
to what they were in the absence of oscillatory external disruption (i.e. for A — 0). 
The effects of impact correction are relatively minor, limited to enlarging the frozen 
phase F, to dampening the fluctuations, and to reducing the non-ergodicity effects in 
the anomalous phase A. 
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7. Discussion 

In this paper we have studied analytically the response of spherical minority games 
(MG) to external time-dependent market disruption of the overall market bid, which 
is envisaged to model the effects of either market regulators or other socio-political or 
natural events outside the market that impact directly on the asset demand and hence 
the asset price. The advantages of the spherical formulation of the game are that one 
can derive fully explicit (and exact) macroscopic dynamical equations and work out the 
market's response rigorously and effectively, even if one includes self-impact correction 
(where agents correct their actions for their own impact on the market), and that the 
volatility can be calculated exactly. We have focused specifically on evaluating and 
comparing the long-time consequences of two types of external bid terms: stationary 
ones, Ae{t) — A, and oscillating ones, Ae{t) — A(— 1)*. These two cases give rise to 
different phase diagrams and different behaviour of the market volatility. The presence 
of an oscillating external bid term is disasterous for the efficiency of the market: the 
market's response is always to oscillate in phase with the perturbation, making the 
situation worse and increasing significantly the volatility. Yet, rather surprisingly, 
the phase diagram in this case is not at all affected by the imposed external bid 
oscillation. All this is more or less the opposite of how the market responds to 
stationary external bid terms, where the phase diagram is strongly affected by the 
amplitude of the bid perturbation, but the fluctuations remain largely unaffected. 
All theoretical predictions have been tested against numerical simulations, and found 
to be supported perfectly well. Our present study emphasizes the need for further 
development of those versions of the generating functional analysis formalism that 
include the stochastic overall bid dynamics explicitly, as in [Til IE] , so that one can 
also tackle the MG market's response to events in its environment for non-spherical 
(i.e. more realistic) versions, where the route to explicit macroscopic equations that 
proved fruitful in this paper is no longer available. 
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